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^ . Abstract 

Loop is a powerful program construct in classical computation, but its power is still not 
exploited fully in quantum computation. The exploitation of such power definitely requires a 
■ deep understanding of the mechanism of quantum loop programs. In this paper, we introduce 

| a general scheme of quantum loops and describe its computational process. The notions of 

termination and almost termination are proposed for quantum loops, and the function computed 
by a quantum loop is defined. To show their expressive power, quantum loops are applied in 
describing quantum walks. Necessary and sufficient conditions for termination and almost 
termination of a general quantum loop on any mixed input state are presented. A quantum 
\Q • loop is said to be (almost) terminating if it (almost) terminates on any input state. We show 

that a quantum loop is almost terminating if and only if it is uniformly almost terminating. 
It is observed that a small disturbance either on the unitary transformation in the loop body 
i[ or on the measurement in the loop guard can make any quantum loop (almost) terminating. 

Moreover, a representation of the function computed by a quantum loop is given in terms of 
finite summations of matrices. To illustrate the notions and results obtained in this paper, two 
simplest classes of quantum loop programs, one qubit quantum loops, and two qubit quantum 
qh loops defined by controlled gates, are carefully examined. 
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1 Introduction 

One of the most striking advances in quantum computing was made by Shor [26] in 1994. By 
exploring the power of quantum parallelism, he discovered a polynomial-time algorithm on quantum 
computers for prime factorization of which the best known algorithm on classical computers is 
exponential. In 1996, Grover [15] offered another apt killer of quantum computation, and he found 
a quantum algorithm for searching a single item in an unsorted database in square root of the 
time it would take on a classical computer. Since both prime factorization and database search 
are central problems in computer science and the quantum algorithms for them are highly faster 
than the classical ones, Shor and Grover's discoveries indicated that quantum computation offers 
a way to accomplish certain computational tasks much more efficiently than classical computation 
and thus stimulated an intensive investigation on quantum computation. After that, quantum 
computation has been an extremely exciting and rapidly growing field of research. In particular, 
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a substantial effort has been made to find new quantum algorithms and to exploit the techniques 
needed in building functional quantum computers. 

Currently, quantum algorithms are expressed mainly in the very low level of quantum circuits. 
In the history of classical computation, however, it was realized long time ago that programming 
languages provide a technique which allows us to think about a problem intended to solve in 
a high-level, conceptual way, rather than the details of implementation. Recently, in order to 
offer a similar technique in quantum computation, some authors begun to study the design and 
semantics of quantum programming languages. In the pool of imperative languages, the earliest 
proposal for quantum programming language was made by Knill in |18j . where a set of basic 
principles for writing quantum pseudo-code was outlined and an imperative pseudo-code suitable 
for implementation on a quantum random access machine was defined. The first real quantum 
programming language, QCL, was proposed and a simulator for this language was implemented 
by Omer \21\ I22j. A quantum programming language in the style of Dijkstra's guarded-command 
language, qGCL, is designed by Sanders and Zuliani in [23j EZl EHl [29] . A probabilistic predicate 
transformer semantics of qGCL was given, a refinement calculus for it was introduced, and a 
compiler from qGCL to a simple quantum architecture was defined. A quantum extension of C+- 1- 
was also proposed by Bettelli et al [3J, and it was implemented in the form of a C++ library. In 
the functional programming style, the first quantum language, QFC, was defined by Selinger [21] 
based on the idea of classical control and quantum data. Programs in the language QFC are 
represented via a functional version of flow charts, and QFC has a denotational semantics in 
terms of complete partial orders of super-operators. In addition, quantum process calculus CQP 
(Communicating Quantum Processes) was introduced by Gay and Nagarajan [HI [14], and QPAlg 
(Quantum Process Algebra) was proposed by Jorrand and Lalire [HI [19] in order to support the 
formal specification and verification of quantum cryptographic protocols. Also, Feng et al [H] 
defined a model qCCS of quantum processes, which is a natural quantum extension of classical 
value-passing CCS with the input and output of quantum states, and unitary transformations 
and measurements on quantum systems. Semantic techniques for quantum computation have also 
been investigated in some abstract, language-independent ways. For example, a notion of quantum 
weakest precondition is introduced and a Stone-type duality between the state transition semantics 
and the predicate transformer semantics for quantum programs is established by D'Hondt and 
Panangaden [8], and proof rules for probabilistic programs were generalized by Feng et al [10] 
to purely quantum programs. There are already two excellent surveys on quantum programming 
languages and related researches [25j [12] 

Loop is a powerful program construct in classical computation [9]. In the area of quantum com- 
putation, looping technique has also attracted a few authors' attention. For example, Bernstein 
and Vazirani [5j [6] introduced some programming primitives including looping in the context of 
quantum Turing machines; some high-level control features such as loop and recursion are provided 
in Selinger's functional quantum programming language QFC. However, the full power of quantum 
loop programs is still to be exploited. The exploitation of such power definitely requires a deep 
understanding of the mechanism of quantum loops. The purpose of this paper is to examine thor- 
oughly mechanism of quantum loops in a language-independent way, and to give some convenient 
criteria for deciding termination of a general quantum loop on a given input. 

This paper is organized as follows. Section 2 is a preliminary section in which some basic 
notions from quantum mechanics needed in this paper are reviewed. In Section 3, a general 
scheme of quantum loop programs is introduced, the computational process of a quantum loop 
is described, and the essential difference between quantum loops and classical loops is analyzed. 
In addition, we introduce the notions of termination and almost termination of a quantum loop. 
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The function computed by a quantum loop is also defined. Quantum walks are considered to show 
the expressive power of quantum loops. In Section 4, we find a necessary and sufficient condition 
under which a quantum loop program terminates on a given mixed input state (Theorem 14. ip . 
In Section 5, a similar condition is given for almost termination (Theorem 15 . X [) . Furthermore, we 
prove that a quantum loop is almost terminating if and only if it is uniformly almost terminating 
(Theorem 15. 2p . and a small disturbance either on the unitary transformation in the loop body 
(Theorem 15 .3|) or on the measurement in the loop guard (Theorem I5.4|) can make any quantum 
loop (almost) terminating. In Section 6, a representation of the function computed by a quantum 
loop is presented in terms of finite summations of complex matrices (Theorem I6.2|) . To illustrate 
the notions and results presented in the previous sections, Sections 7 is devoted to some examples 
which observe the computational behavior of two simplest classes of quantum loops: one qubit 
loops, and two qubit loops defined by controlled operations. Section 8 is the concluding section in 
which we draw the conclusion and point out some problems for further studies. 



2 Preliminaries 

For convenience of the reader we briefly recall some basic notions from quantum theory and fix the 
notations needed in the sequel. We refer to [20] for more details. 

An isolated physical system is associated with a Hilbert space which is called the state space 
of the system. We only need to consider finite dimensional Hilbert space in quantum computation. 
An n— dimensional Hilbert space is an n— dimensional complex vector space H together with an 
inner product which is a mapping (-|-) : H X H — ► C satisfying the following properties: 

1- > with equality if and only if |0) = 0; 

2. (0|^) = 

3. (0|AiV>i + \21P2) = Ai(0|Vi> + A 2 (0|^ 2 >, 

where C is the set of complex numbers, and A* stands for the conjugate of A for each complex 
number A E C. For any vector in H, its length \\\tp)\\ is defined to be yj Let V be a set 
of vectors in a Hilbert space H. Then span{V) is defined to be the subspace of H spanned by V, 
that is, it consists of all linear combinations of vectors in V. An orthonormal basis of a Hilbert 

{1, if i = j, 
Then the trace of a linear operator A on H 
0, otherwise. 

is defined to be tr{A) = ^2i(i\A\i). 

A pure state of a quantum system is a unit vector in its state space, that is, a vector \tp) with 
|||-0)|| = 1, and a mixed state is represented by a density operator. A density operator in a Hilbert 
space H is a linear operator p on it fulfilling the following conditions: 

1. p is positive in the sense that (ip\p\ip) > for all \if>); 

2. tr(p) = 1. 

An equivalent concept of density operator is ensemble of pure states. An ensemble is a set of 
the form {(pi, \ipi})} such that Pi > and \ipi) is a pure state for each i, and YliPi = 1- Then 
p = 'Yl i Pi\ipi}{ipi\ is a density operator, and conversely each density operator can be generated by 
an ensemble of pure states in this way. The set of density operators on H is denoted T>(H). A 
positive operator A is called a partial density operator if tr(A) < 1. We write T>~(H) for the set 
of partial density operators on H. Obviously, T>(H) C V~{H). 
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The evolution of a closed quantum system is described by a unitary operator on its state space. 
A linear operator U on a Hilbert space H is said to be unitary if U'U = Ih, where Ih is the 
identity operator on H, and W is the adjoint of U. If the states of the system at times ti and i 2 
are p\ and P2, respectively, then p2 = U p\U^ for some unitary operator U which depends only on t% 
and ti- In particular, if p\ and P2 are pure states \ipi) and respectively, that is, pi = 
and P2 = 1^2) {ifa], then we have \ip2) = U\ipi). 

A quantum measurement is described by a collection {M m } of measurement operators, where 
the indexes m refer to the measurement outcomes. It is required that the measurement operators 
satisfy the completeness equation Y^m Mrn,M m = Ih- If the system is in state p, then the probability 
that measurement result m occurs is given by p(m) = tr(MmM m p), and the state of the system 
after the measurement is M p^ m ■ For the case that p is a pure state that is, p = we 
have p(m) = ||M m |^)|| 2 , and the post-measurement state is -^==. In particular, a projective 

measurement is described by an observable which is represented by a Hermitian operator. A 
Hermitian operator is a linear operator M with = M. An eigenvector of a linear operator A 
is a nonzero vector |A) such that A\X) = A | A) for some A £ C, where A is called the eigenvalue 
of A corresponding to |A). We write spec(A) for the set of eigenvalues of A which is called the 
spectrum of A. It is well known that all eigenvalues of a Hermitian operator M are reals. Let 
M = Ylrnespec(M) ra ^m be the spectral decomposition of M where for each m £ spec(M), P m 
is the projector to its corresponding eigenspace. Obviously, these projectors form a quantum 
measurement {P m : m £ spec(M)}. If the state of a quantum system is p, then the probability that 
result m occurs when measuring M on the system is p(m) = tr(P m p), and the post-measurement 
state of the system is F ™^™ ■ 

The state space of a composite system is the tensor product of the state spaces of its components. 
Let Hi and H2 be two Hilbert spaces. Then their tensor product Hi <g> H2 consists of linear 
combinations of vectors \1p1ip2) = l^i) ® 1*02) with \ipi) £ Hi and \ip2) £ H2. For any linear 
operator Ai on Hi and A2 on H2, Ai ® A2 is an operator on Hi ® H2 and it is defined by 

(A x ® A 2 )|^i^2> = ® Vl2|^2> 

for each £ -ff 1 and ^2) £ ^2- Since density operators are special linear operators, their tensor 
product is then well-defined. If component system i is in state pi for each i, then the state of the 
composite system is (^iPi- 



3 Basic Definitions 

We first give a general and formal formulation of quantum loop programs. Suppose that we have n 
quantum registers qi, . . . , q n , and their state spaces are Hi, ... , H n , respectively. We further assume 
that U is a unitary operator on the tensor product space H = ®" =1 Hi. Let M = ^2 m mP m be a 
projective measurement on H. Then for any X C spec(M), the quantum loop program defined by 
U, M and X may be written as follows: 

while (M £ X) {q := Uq} (1) 

where q is used to denote the sequence qi,...,q n of quantum registers. Let Px = Ylm&x^m- 
and P^ = Ih — Px = Yl m espec(M)-X ^rm where Ih is the unit operator on H. Then the guard 
"M £ A" of loop ([T]) means that the projective measurement {Px,Px~} 1S applied to q, and the 
outcome corresponding to Px is observed. The body of the loop is the assignment il q := Uq" , that 
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Figure 1: Quantum loop ([I]) 



is, a command of performing unitary transformation U on the state of the sequence q of quantum 
registers. This loop can be visualized by Fig.l. 

It is worth noting that the projective measurement we perform to check the guard condition of 
loop ([I]) is {Px, P~x} rather than M itself, because we need only tell whether or not the measurement 
outcome belongs to X. Any further information about the exact outcome is useless, and will bring 
unnecessary disturbance to the system we measured. 

We now examine the computational process of the above loop program. For any input state 
Po = p £ T>(H), the behavior of the above quantum loop can be described in the following unwound 
way (see Fig. 2): 

1. This is the initial step. The loop program performs the projective measurement {Px,Px~} 
on the input state p. If the outcome corresponding to Px is observed, then the program per- 
forms the given unitary operation U on the post-measurement state. Otherwise the program 
terminates. Formally, the loop will terminate with probability Pj-\p) = t r {P~x~P) anc ^ ^ wm 
continue with probability p^L(p) = 1 — p^\p) = tr{Pxp)- In the case of termination, the 
output of the loop is 

fi) = p xP p x 
P ° Ut 

and in the case of nontermination, the state of q system after the measurement is 

(i) _ PxpPx 
^ " &(p) ' 

Furthermore, p^ id will be fed to the unitary operation U and then the state p^ = U p^ id U^ 
is returned, which will be used as the input state in the next step. 

2. This is the induction step. Suppose that the loop has performed n steps, and it did not 

(n) (n) _ 

terminate at the nth step, that is, p N ^ > 0. If p\ n is the state of q system returned at the 

nth step, then in the (n + l)th step, the termination probability is Pt (p) = ^i^xPin) 
and the output is 

p n i n )p 
(n+l) = f_x£in fJL 

P ° Ut p^ +l \p) ' 

The loop continues to perform the unitary operation U on the post-measurement state 

(n+l) _ PxP^Px 

Pmid 



Pnt 1] (p) 
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Figure 2: The computational process of loop ([!]) 



with probability ' (p) = 1 - p^ '(p) = tr(P x pln), and the state = Up^ 'U^ 

will be returned. 

Note that not only a pure quantum state but also a mixed state is allowed to feed into a 
quantum loop. In fact, quantum programming with mixed states has already been considered in 
the previous literature; for example, see [MIES]- 

There is an essential difference between the computing process of quantum loops and that of 
classical loops. In a classical loop the states of variables do not change during verification of the 
loop condition. However, in a quantum loop it is impossible to check the loop condition directly. 
Instead, the loop program needs to extract information about the registers q%, . . . , q n by performing 
a measurement M and thus their states will be changed. 

To demonstrate the expressive power of quantum loops, let us consider an interesting example. 
Quantum walk is a natural quantum extension of classical random walk, which in turn has proved 
to be a fundamental tool in computer science, especially in the designing of algorithms [17]. In this 
example, we consider a discrete coined quantum walk on an n-cycle with an absorbing boundary 
at position 1, and express this kind of quantum walk in our language of quantum loops. For 
more details about quantum walk on a cycle, or more generally, on any graph, we refer to pQ. 
The following example shows that a quantum walk can be described very well in the language of 
quantum loops. 

Example 3.1 Let Ha be a 2- dimensional 'coin' space with orthonomal basis states |0) and |1), 
and Hy be the n- dimensional principle space spanned by the position vectors \i) : i = 0, . . . ,n — 1. 
Then each step of the quantum walk we are concerned with consists of three sub-steps: 

1. A 'coin-tossing operator' H = |+)(0| + |— }(1| is applied to the coin space, where |+) = 
(|0) + |1»/V2 and |-) = (|0) - |l»/\/2. 

2. A shift operator 

re— 1 Ti—l 

fir = ^ \i G ® |0)(0| + \i © |1)(1| 

i=0 i=0 

is performed on the space Hy <g> Ha, which makes the quantum walk one step left or right 
according to the coin state. Here and © denote subtraction and addition modulo n, respec- 
tively. 

3. Measure the principle system to see if the current position of the walk is 1. If the answer is 
'yes' then terminate the walk, otherwise the walk continues. 

Formally, we can formulate the walk described above by a quantum loop: 

while (M / 1) {q := Uq} (2) 

where M = Y27=o ® I2, U = S{I n (8> H), and q is a quantum register in Hy Ha- For 

simplicity, we write M 7^ 1 in the loop guard to denote M £ X for X = {0, 2, . . . , n — 1}. 

One of the most important problems concerning the behavior of a loop program is its termina- 
tion. 

Definition 3.1 1. If p^ T {p) = for some positive integer n, then it is said that the loop with 
input p terminates. 
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2. The nonterminating probability of the loop with input p is defined to be 

Pnt(p) = lim PntQ 3 ) 

where (and in the sequel) p 1 ^ T (p) = Y17=xPnt(p) denotes the probability that the loop does 
not terminate after n steps. 

3. We say that the loop with input p almost terminate whenever Pnt(p) = 0. 
4- If Pnt(p) > 0, then we say that the loop with input p does not terminate. 

Intuitively, a quantum loop almost terminates if for any e > 0, there exists a big enough positive 
integer n(e) such that the probability that the loop terminates at the n(e)th step is greater than 
1 — e. Obviously, if a quantum loop terminates on a given input state, then it also almost terminates 
on the same input. 

A classical loop may terminate or not, but a quantum loop has an additional possibility of 
almost termination. Clearly, this is caused by the probabilistic nature of quantum mechanics. 

Definition 3.2 1. A quantum loop program is said to be terminating (resp. almost terminating) 
if it terminates (resp. almost terminates) with all input p € T>(H). 

2. A quantum loop is uniformly almost terminating if for any e > there exists a positive integer 
n(e) such that p r ^ T (p) < e holds for all input p € T>(H) whenever n > n(e). 

It is clear that uniformly almost terminating quantum loops are almost terminating. 
Note that the case of X = or spec(M) is trivial. In fact, the loop (pQ) is equivalent to 

while (false) {q := Uq} 

when X = 0, and it is equivalent to 

while (true) {q := Uq} 

when X = spec(M). The former terminates immediately and does nothing, and the latter will 
loop forever. In what follows we always assume that C X C spec(M). 

In the computational process of a loop program, a density operator is input, and a density 
operator is outputted with a certain probability at each step. Thus, we have to synthesize these 
density operators returned at all steps according to the respective probabilities into a single one as 
the overall output. Note that sometimes the loop does not terminate with a nonzero probability. 
The synthesized output may not be a density operator but only a partial density operator. Then 
a loop defines a function from density operators to partial density operators on H. 

Definition 3.3 The loop |7p defines a function F : V(H) — > V~(H) in the following way: 

oo 

F(p) = Pout = ^PNT 1)+ (p)PT } (p)Pout 
n=l 

for each p E D(H). The function F is called the function computed by the loop {2JJ. 
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It should be noted that in the defining equation of F(p) the quantity P^t^ ~(p)Pt'\p) ^ s the 
probability that the loop does not terminate at steps from 1 to n — 1 but it terminates at the nth 
step. 

For the case that p is a pure state, that is, p = \ip){tp\ for some pure state \ip), we will write 
F(\tp)) in place of F(p) for simplicity. 

In the remainder of this section, we are going to present some basic properties of quantum 
loops. For any operator A on H, we write Ax = PxAPx, that is, Ax is the restriction of A on the 
subspace of H corresponding to the projector Px- First, the computational process of quantum 
loop (pQ) can be summarized as: 



Lemma 3.1 Let p be the input state to the loop UP- Then for any positive integer n, we have: 



p n N + T (p)=tr(U x -'p x U^) 



n-l 



An-U 



and 



F(p) = P x pP x + P T U [J2 U xPxU' x n WP X . 

\n=0 / 



Proof. First, it is easy to check by induction on n that 



Pnt(p) 



Hu^pxU^- 1 ) 

tr(U x - 2 P xU^ 2 ) 
tr(px), 



if n > 2, 



if n = 1, 



(3) 
(4) 



(5) 



and 



P { t\p) 



tr{P x UU x - 2 Px UT 2 U X ) 
tr(U x - 2 p x U^- 2 ) 



1 - tr(px), 



if n > 2, 



if n = 1, 



.(») 

Pout 



P x UU x - 2 p x U^- 2 WP x 



tr(P x UU x - 2 p x Up~ 2 W) 



(6) 



Then Eq. ([3]) follows from Eq. ([5]), and Eq. (@J comes from Eqs. ©, ©, and ([7|). 



(7) 
□ 



It should be pointed out here that convergence of the infinite series in Definition 13.31 and Eq. ^) 
is guaranteed by the facts that the set T>~{H) is a directed complete poset under the Lowner order 

{k t 1 00 

^2 n= oU x pxUjP > is non-decreasing in this order. For the details, we refer 

to [23]. On the other hand, from Eq. Q) and Kraus representation theorem (|20j. Theorem 8.1) 
we notice that the function F computed by loop ([T]) is a super-operator (also called quantum 
operation). 

Let Hx be the subspace of H with projector Px, and the subspace with projector P^. The 
following proposition clarifies the range of the function F computed by the loop ([T]). 



Proposition 3.1 For each p £ T>(H), we have: 

1. (4>\F(p)\4>) =0 if\<j>) or\4>)eH x ; 
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MX) 



N 



( X N 





V o 



A 



jV-1 



X N 







JV 
2 

N 
1 



X N-2 
X N-1 








N 
r-2 

N 
r — 3 

JV 
r - 4 



^N-r+2 
X N-r+3 



X N 




r - 1 

JV 
r-2 

JV 
r - 3 



JV 



\AT-r+2 
^iV-r+3 

X N-1 



• (8) 



J r (A)^v = (£(") ^vi+i, E ( T ) ^+2, - • ^"r-i + ( T ) ^ V ' A ^) T - (9) 
i=o ^ 1 ' i=o ^ 1 s V / 



<?. tr(F(p)) = 1 — pnt(p)- Thus, F(p) G T>(H) if and only if the loop |7p uflft/i mpwi staie p 
almost terminates. 

Proof. 1. By definition we know that fM<?!>) = ^M^) = ^- Then it follows immediately from 
Lemma 13.11 

2. By induction on A; it may be shown that 

&fe 1)+ 0»)#W = 1 

n=l 

Then we have: 

oo 

«=Epfe 1)+ (c)!4 fl) w 

ra=l 

= 1 - lim f$r(» = 1 - Pnt(p)- 

n— >oo 

The conclusion follows immediately. □ 
From Proposition 13.11 1 we see that F is indeed a function from T>(H) into D~{H-y). 

4 Termination 

The aim of this section is to give a necessary and sufficient condition under which the loop (JT]) 
terminates on a given input state. 

We first give a lemma which allows us to decompose an input density matrix into a sequence 
of simpler input density matrices when examining termination of a quantum loop. 

Lemma 4.1 Let p = YliPiPi where Pi > and p,i S T>(H) for each i, and J2i.Pi = 1- Then the 
loop HP with input p terminates if and only if it terminates with input pi for all i. 
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Proof. For each i, if the loop (JTJ) with input pi terminates, then there exists a positive integer rjj 
such that p^^(pi) = 0. Let no = maxj n^. Then p^^(pi) = for all i, and this yields 

i 

Conversely, if the loop (pQ) with input p terminates, then there exists a positive integer no such 
that p^^{p) = 0. This implies that for each i, p^^{pi) = because p^T(pi) > for each i. □ 

If {(pi, \ipi))} is an ensemble with pi > for all i, and p = ^2iPi\ipi){ipi\, then the above lemma 
asserts that the loop (pQ) terminates on input mixed state p if and only if it terminates on input 
pure state for all i. In particular, we have: 

Corollary 4.1 A quantum loop is terminating if and only if it terminates with all pure input 
states. 

Second, the termination problem of a quantum loop may be reduced to a corresponding problem 
of a classical loop in the field of complex numbers. Let |mi), |ni2), • • • , \itik) be an orthonormal 
basis of H such that 

k K 

yjm^jmi] = P x , y~] \mi){mi\ = Px, 

i=l i=fc+l 

where 1 < k < K. Without any loss of generality, we assume in the sequel that the matrix repre- 
sentations of U, Ux,px (denoted also by U, Ux,px respectively for simplicity) are taken according 
to this basis. In particular, for each pure state we write \ip)x f° r the vector representation of 
Px\ip) under this basis. 

Lemma 4.2 The quantum loop {!]) terminates with input p £ D(H) if and only if PxUj^ = 
Ofcxfc f or some nonnegative integer N , where Ofe X fc is the (k x k)-zero matrix. In particular, it 
terminates with pure input state if and only if U^\ijj)x = for some nonnegative integer N, 
where is the zero vector of length k. 

Proof. This result follows from Eq. ([3]) and the fact that tr{A) = if and only if A = when A is 
positive semi-definite. □ 

Third, we show certain invariance of termination of a classical loop under a nonsingular trans- 
formation. 

Lemma 4.3 Let S be a nonsingular (k x k) — complex matrix. Then the (classical) loop: 

while (v / 0) {v := U x v} (v G C k ) 
terminates on input vo £ C k if and only if the following loop: 

while (v + 0) {v := {SUxS'^v} (v £ C k ) 

terminates on input Svq. 

Proof. Note that SV ^ if and only if v ^ because S is nonsingular. Then the conclusion follows 
from a simple calculation. □ 

Furthermore, we shall need the famous Jordan normal form theorem in the proof of the main 
result in this section. 
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Lemma 4.4 (Jordan normal form; ^) For any (k x k) — complex matrix A, there is a nonsingular 
(k x k) — complex matrix S such that A = SJ(A)S~ 1 , where 



J (A) = diag(J kl {\i), Jfc 2 (A 2 ), . . . , <4,(Az)) 

( ^(Ai) \ 
^2(^2) 



V 




\ 



(10) 



1 



V 



A; / 



is a {ki x ki)- Jordan block for each 1 < i < I. Furthermore, if the Jordan blocks corresponding to 
each distinct eigenvalue are presented in decreasing order of the block size, then the Jordan normal 
form is uniquely determined once the ordering of the eigenvalues is given. 

The following technical lemma is also needed. 

Lemma 4.5 Let J r (A) be a (r x r) — Jordan block, and v G C r . Then J r (\) N v = for some 
nonnegative integer N if and only if A = or v = 0, where is the zero vector of length r. 

Proof. The "if part is clear. We now prove the "only if part. By a routine calculation we obtain 
the matrix J r (X) N as in Eq. ([8]). Notice that J r (X) N is an upper triangular matrix with the diagonal 
entries being A^. So if A 7^ then J r (X) N is nonsingular, and then J r (\) N \r = implies v = 0. □ 

Now we are able to present the main result of this section. 
Theorem 4.1 Suppose the Jordan decomposition of Ux is Ux = SJ(Ux)S~ 1 , where 



Let S' _1 |V')x be divided into I sub-vectors vx,V2, . . . ,v; such that the length of Vj is fcj. Then the 
loop (Jl) terminates on input if and only if for each 1 < i < I, Aj = or Vj = ; where is the 
zero vector of length ki . 

Proof. Using Lemmas 14.21 and 14.31 we know that the loop ([T|) terminates on input if and only 
if J(Ux) N S~ 1 \t/j)x = for some nonnegative integer N. A simple calculation yields 



Therefore, J(Ux) N S 1 \ip)x = for some nonnegative integer N if and only if for each 1 < i < I, 
there exists a nonnegative integer Ni such that Jki{\) Ni ^i = 0. Then we complete the proof by 



J(Ux) = diag(J kl (\i), Jfe 2 (A 2 ), ... , Jfej(Aj)). 



JiUx^S-^x 



((JkM N Vl) T , (Jta(A 2 )*V 2 ) T , . . . , (J km (\m) N V m ) T ) T . 



using Lemma 14.51 



□ 
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Corollary 4.2 Loop ^ is terminating if and only if Ux has only zero eigenvalues. 

Theorem 14.11 gives a necessary and sufficient condition for termination of loop (pQ) on an input 
pure state. Obviously, we can decide whether the loop ([T]) terminates on any given mixed state as 
input by combining Lemma 14.11 and Theorem 14.11 The condition for termination of loop ([I]) can 
be considerably simplified in the special case when Ux is normal, that is, UxU x = U x Ux- In this 
case, Ux has the following simple spectrum decomposition: 

k 
i=l 

Then from Eq. (|3|) we have for any p £ D(H): 

k 

pVrip) = me ^wmmr 1 ) 

i=l 

This implies immediately the following: 

Corollary 4.3 Suppose Ux is normal and its spectrum decomposition is given by Eq. Ul\) . Then 
we have: 

1. loop (QP terminates on input state p if and only if for any i = l,...,k, Aj ^ implies 
{i\p\i) = 0, or equivalently, tr(Uxp) = 0. 

2. loop ([T]j is terminating if and only if Ux = 0. 



(11) 



(12) 
(13) 



5 Almost termination 

In this section we are going to present a necessary and sufficient condition under which the loop JT]) 
almost terminates on any given input state. We first give a lemma similar to Lemma |4. II so that a 
mixed input state can be reduced to a family of pure input states. 

Lemma 5.1 Let p = 'YlnPiPi where Pi > and pi E T3(H) for each i, and J2i.Pi = 1- Then the 
loop HP with input p almost terminates if and only if it almost terminates with input pi for all i. 

Proof. Notice that Px^pip) = Yli PiP^rr(Pi) from Eq. ([3]). The result follows by taking limits about 
n in both sides of the above equation. □ 

Corollary 5.1 A quantum loop is almost terminating if and only if it almost terminates on all 
pure input states. 

The following lemma is a key step in the proof of our main result in this section. 

Lemma 5.2 The loop {2P almost terminates on the pure input state \ip) if and only i/linin^oo ||C/^|?/')|| = 
0. 
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Proof. PromEq. © wehavep^! r (|V')) = \\U X 
0. 



n-li 



! . Soptvt(IV')) = if and only lim r , 



\U^)\\ 
□ 



The following theorem gives a necessary and sufficient condition for almost termination of a 
quantum loop on a pure input state. 



Theorem 5.1 Suppose that Ux, S, J{Ux), Jki(^i) andvi (1 < i < I) are given as in Theorem \4-l\ 
Then the loop (OP almost terminates on input \tp) if and only if for each 1 < i < I, |Aj| < 1 or 
Vj = 0, where is the zero vector of length k{. 



Proof. First, for any nonnegative integer n, we have U x 
if and only if 

lim WJiUxTS' 1 



SJiUxTS- 1 ] 
= 



Then lim r 



\jjn 



(14) 



since S is nonsingular. Using Lemma 15.21 we know that the loop (pQ) almost terminates if and only 
if Eq. (HU) holds. Note that 

JiUxTs- 1 ^) 

= ((J kl (\i) n vi) T , (J k2 (A 2 rv 2 ) T , . . . , ( J kt (A,)"vO T ) T - 



Then Eq. (I14j) holds if and only if 



lim || J ki (\i)vi\ \ = 

n— >oo 

for all 1 < % < I. Furthermore, for each 1 < i < I, from Eq. 
only if the following k{ equations are valid: 



(15) 

we see that Eq. (|15p holds if and 



lim r , 
lim r , 



fcj— 1 


1 n 


j=Q 


\ J 




ki-2 


( n 


j=0 


\ 3 





\ n ~l n 



\r 3 v. 



lim,. 



Pn v i{k-l) + 



n 
1 



(16) 



hm. 



A>i fc = 0, 



where it is assumed that Vj 
If |AJ < 1, then linin^oo 



(va,Vi2, ■ ■ .,v iki ). 



n 



J 



X 



n-j 



for any < j < ki — 1, and all of the above equations 



in Eq. (|16p follow. On the other hand, if |Aj| > 1, then from the last equation in Eq. (|16p we 
know that v^. = 0. Putting Vi k = into the second equation from bottom in Eq. (|16|) we obtain 
v i(k— l) = 0- We can further move from bottom to top in Eq (|16p in this way, and finally we get 
v-ii = Vi2 = • • • = = = 0. This completes the proof. □ 

Corollary 5.2 Loop (OP is almost terminating if and only if all the eigenvalues of Ux have norms 
less than 1. 

In the case when Ux is normal, we have the following corollary which is also easy to prove 
directly from Eq. (|13|) . 
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Corollary 5.3 Suppose Ux is normal and Eq. ill]) is its spectrum decomposition. Then 

1. loop (QP with input p almost terminates if and only if for any i = 1, . . . , k, |Aj| = 1 implies 
{i\p\i) = 0, i.e., the set 

I' = {l<i<k \ \Xi\ = 1 and (i\p\i) > 0} 

is empty. 

2. The nonterminating probability is Pnt(p) = J2i£i'(i\p\i) • 

Now we are able to show that the quantum walk ([2]) in Example 13.11 is almost terminating. It 



is 



direct to calculate that Px = !*)(*! ® I2 an d 



u x = Yl (NX*© 1 ! ® +|*©1><*| ®|0)<-|)- 

With Corollary 15.21 it suffices to prove that each eigenvector of Ux has its norm strictly less that 1. 
By contradiction, suppose Ux has an eigenvalue A with unit norm, and one of the corresponding 
normalized eigenvector is 

n-l 



|V) ='52(a i \i)®\+)+Pi\i)®\-)), (17) 

i=0 

where £i(K| 2 + \Pi\ 2 ) = 1. Then we have 



A|V) = U x \i/>) = ( a «wl*> ® 1°) + Al* © : > ® I 1 ))- ( 18 ) 



Comparing Eqs. (|17p and (|18p . we derive further that 



and 



J °« if i = 0,l; , . 

[ V2a ie i/A, if 1 + 0, 1 

a ^ A = i v^AexA, if Ml, 2. (20) 
On the other hand, since |A| = 1, we know from Eq. (118p that 

^(|a iel | 2 + |A| 2 ) = ||A|V)|| 2 = l. (21) 

So we have: 

ax = a 2 = Po = A = 0. (22) 

Taking Eq. ()22[) back into Eqs. (|19p and (|20p we can deduce that ai = Pi = for any i. This is a 
contradiction. 

To conclude this section, we observe some further properties of almost terminating quantum 
loops. The following theorem indicates that the notion of uniformly almost terminating loop 
coincides with almost terminating loop. 

Theorem 5.2 The quantum loop ([I]) is almost terminating if and only if it is uniformly almost 
terminating. 
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Proof. If loop dU) is almost terminating, then we have |Aj| < 1 for any i = 1, . . . , I from Corollary 
15.21 Let Ux = SJ(Ux)S~ 1 be the Jordan decomposition of Ux- Then from Eq. ([3]) we have: 

p n N + T (p) = \\SJ(Ux) n - 1 S~ 1 pf\\ 2 . 

By using the properties of matrix norm, we derive that 

PnMp) < (II^MI^II-llpf \\) 2 \\J(U x )f^ 

< (iisii-ii^- 1 !!) 2 !!^)!! 2 ^- 1 ). 

Since spec(J(Ux)) G [0, 1), from Eq. (jHJ) we can check easily that || J(?7x)|| n — ► when n — ► oo. So 
for any e > 0, we can take n(e) large enough such that 



\\J(Ux)\\ 



2(n(e)-l) 



< 



15" 



-Hh2 



Then we have p r N ^p(p) < e for all p whenever n > n(e). Thus loop (pQ) is uniformly almost 
terminating. □ 

The next two theorems show that the notion of almost terminating loop is sensitive. More 
explicitly, a small disturbance either on the unitary transformation in the loop body or on the 
measurement in the loop guard can make any quantum loop almost terminating. 

We first need to prove a technical lemma. 



Lemma 5.3 Suppose 
norm. Then: 



\i) is an eigenvector of Ux and its corresponding eigenvalue Aj has unit 



1. \i) G Hx, and it is also an eigenvector of U with an eigenvalue of unit norm; 

2. P Y U\i) = 0. 

Proof. Assume that Ux\i) = A|z) and |A| = 1. First, we see that \Px\i) = PxUx\i) = Ux\i) = A|i). 
Thus, P X \i) = \i) and \i) G Hx- Furthermore, by the Gram-Schmidt procedure we can find an 
orthonormal basis {\j}} for H which contains We assume that U\i) = ^ • Pj\j) and Ylj \l x j\ 2 = 
1. Then it holds that \m\ = \{i\U\i)\ = \(i\P X UP x \i)\ = |A| = 1. This implies that pj = for all 



j ^ i, and U\i) = pi\i). Finally, PxU\i) = p { P- 



x 



0. 



□ 



Theorem 5.3 For any M , X ^ spec(M) and U in loop HP, and for any e > 0, there exists a 
unitary operator V such that \\U — U'\\ < e and the following loop is almost terminating: 

while (M G X) {q := U'q}. 

Proof. By using Corollary 15.21 we only need to find a unitary operator U' such that ||J7 — U'\\ < e 
and all eigenvalues of PxU'Px have norms less than 1. On the other hand, Lemma 15.31 implies 
that a necessary condition for PxU'Px to have an eigenvalue with unit norm is that U' has an 
eigenvector lying in the space Hx- Here Hx is the subspace with projector Px- So we need 
only to show that we can take U' such that ||J7 — U'\\ < e and at the same time none of the 
eigenvectors of U' lies in Hx- To achieve this, we first write out the spectrum decomposition of U 
as U = Y^i f^i\'4 , i)('ipi\- If each Hx then we have done. Otherwise suppose \ipi ) G Hx for 
some io- From X / spec(M) there exists jo such that \ipj ) ^ Hx- Let 



(23) 
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14,} = \/T^IV>; > - rfWio), (24) 

and |i/^) = for i ^ ioijo- Here 5 is a very small but positive real number which will be 
determined later. It is obvious that the set \ip'j) are also orthonormal, and \ip'i ), Wj ) Hx- 
Let Ui = Yli A t «IV'i)(V'il- Then the number of eigenvectors of U\ which lie in Hx is strictly less 
than that of U . Repeating the above steps we can finally find a sequence of unitary matrices 
U = Uo, Ui, . . . , Ud, d < K = dim(H), such that all the eigenvectors of Ud does not lie in Hx- 
Take U' = Ud and notice that we can take 5 small enough at each step such that \\Ui— < e/K. 
It then follows that \\U - U'\\ < £to Wi - U i+X \\ < e. □ 

Theorem 5.4 For any M, X ^ spec(M) and U in loop £2]j, and for any e > 0, there exists an 
observable M' with spec(M') = spec(M) such that \\M' — M\\ < e and the following loop is almost 
terminating: 

while (Af G X) {q := Uq} 

Proof. Similar to the proof of Theorem 15.31 it suffices to find M' such that spec(M') = spec(M) , 
|| M — M'\\ < e, and none of the eigenvectors of U lies in H' x , where H' x is the eigenspace of M' 
with eigenvalues in X. Let {\rrii)}fL 1 be an orthonormal basis of H such that Px = \ m i}( m i\- 
Since X ^ spec(M), we have 1 < k < K . 

Let U = ^j\tpj){tpj\ be the spectrum decomposition of U. If all \tpj) Hx then we have 
done. Otherwise assume \ipj ) £ Hx- Then there exists iq < k such that {m-i \ipj ) ^ 0. We take 
some k + 1 < i\ < K and put 

,1 



)=VT^d\m io ) + y/S\m h ), (25) 

) = y/T^5\m il )-y/6\m i0 ), (26) 
where it is required that < S < 1, and \m\) = \rrn) for i ^ io,i%. It is easy to check that the 
set {\m\)}f = i is also an orthonormal basis of H. We write P\ = Xa=i \ m i)( m i I- Let H x be 
the subspace of H with projector P x . Then \tj)j Q ) H x because (mj^ipjg) = —vS(mi \iJjj ) ^ 0. 
Furthermore, we can choose 8 carefully such that the eigenvectors of U in H x is strictly less than 
that in Hx- Indeed, if \iftj) Hx but E H x then it must hold that 

= (ml\ip 3 ) = VT^Simi^j) - v^mJV;)- (27) 

Thus, there are only finitely many S which does not meet our requirement. 

Repeating the above steps we can finally find a sequence of orthonormal bases {|m'),i = 
1, . . . ,K}f =0 with d < K such that |m?) = |m,) for each i < K, and all the eigenvectors of U do 
not lie in H x , the subspace of H with projector P x = XL=i \ m i)( rn i\- Let 

P' m = E \mt)(mf\,M'= Y, mP 'm- 

Pm\m,i)=\mi) m£spec(M) 

Notice that we can take 5 small enough at each step such that 

max{||P! - pj+\ ||P-! - 4 +1 ||} < ' 
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It then follows that 

\\M-M'\\ < £|m|- \\P m -P' m \\ 

m 

< EHfEn^-^n) 

m \i=0 / 

/d-1 \ 

< Eh E 



< e. 



□ 



6 The Function Computed by a Quantum Loop 

In this section, we are going to give a representation of the function computed by the loop (pQ). 
First of all, we consider the simple case that Ux is normal. 

Theorem 6.1 Suppose Ux is normal and its spectrum decomposition is given by Eq. Ul\) . Then 
the function computed by loop is as follows: 

Hp) = PxpPx+ E i^W)01^ 

i,jel 1 3 

where I = {i \ 1 < i < k , \Xi\ < 1}. 

Proof. For any n > 0, we have from Eq. (Ilip that 

k 

P-xUU x = J2>H P x U \ i )(i\ =Y. X * P xU\i)(i\ 

i=l iel 

where the second equality is due to Lemma 15.31 2. Taking this equation into Eq. ([3]) we have: 
F(p) =PxpPx+ E (E*?*r) ^\Px\3)PxU\m^Px- 

i,jel \n=0 ) 

Then the result follows by using Lemma 15.31 1. □ 

We now turn to consider the general case where Ux is not necessary to be normal. To this end, 
the following lemmas are needed. 

Lemma 6.1 (Schur's unitary triangularization; [1]) Given (kxk)-complex matrix A with eigenval- 
ues Ai, . . . , Afc in any prescribed order, there exists a (kxk) unitary matrix V such that A = VTV\ 
where T is upper triangular with diagonal entries Ta = A,, i = 1, . . . , k. 

Lemma 6.2 Let Ux = VTV^ be the Schur's triangularization of Ux- Then for any 1 < i < k, if 
\Tu\ = 1 then Tij = Tji = for all j ^ i. 

Proof. To prove this lemma, we need only to notice T^T = V^U^-UxV < I, and so for any i, the 
Euclidean norms of the i-th row and the i-th column of T must be less than or equal to 1. □ 
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Lemma 6.3 For each Jordan block J r {\) in the Jordan normal form ofUx, if\M = 1> then r = 1. 
That is, each Jordan block corresponding to unit norm eigenvalues has size 1. 



Proof. Suppose Ux = VTV' is the Schur's triangularization of Ux, and the diagonal entries of T 
have been arranged in decreasing order of their norms, i.e., 1 = \Tn\ = ■ ■ ■ = |T#| > |TJ 
• • • > l^fcfcl f° r some t. Then from Lemma 16,2} T must have the form 



m,ml > 



\ 



\ 



T ) 



where T' is (k — t) x (k — t)— dimensional and none of its eigenvalues has unit norm. Let T 1 
S' J(T')S'~ l be the Jordan decomposition of T' , and let 



/Til 



S = V 



It 






S' 



V 



\ 



J{T>) ) 



It is easy to check that SJS 1 is the Jordan decomposition of Ux- Then the result holds from the 
uniqueness of Jordan normal form in the sense presented in Lemma l4.4i □ 

Lemma 6.4 Let J r (A) be a (r x r) — Jordan block, |A| < 1, and v = [v\,V2, ■ ■ ■ ,v r ) T S C r . Then 

oo r— 1 t— 2 



5>(A)»v = £>(AM +1 ,J>(A) 



Vi+2, 



n=0 



i=0 



i=0 



(28) 



/o(A)iy_i + fi(X)v r , fo(X)v r ) 



T 



where 



fi(x) 



d\l - X 
i\ dx 1 



-l 



Proof. For any 1 < m < r, we can see from Eqs. © and ([9]) that the m— component of vector 
E^oWv is 



• uxj \ r—m oo / \ 

Ewv = EE I 

\n=0 / m i=0 n=0 V 7 



,J i+m 



i=0 71=0 

1 — m / ,j ^oo n 



E 

i=0 
r—m 

8=0 



^E^=o^ 



i+m • 



x=A 



The convergence of the above series is guaranteed by the assumption that |A| < 1. 
Now we are able to present the main result of this section. 



□ 
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Theorem 6.2 Suppose that S, J(Ux), ^(Aj) an d Vi(l<i<0 are given as in Theorem \4-l\ 
Without loss of generality, we assume that the Jordan blocks of J(U X ) have been arranged in the 
decreasing order of |Aj|, i.e. 1 = |Ai| = • • • = \Xt\ > |At+i| > • • • > |A/|. Then the output F(\ip)) of 
the loop (OP with input is a (K — k)- dimensional vector lying in the subspace Hx~ : 

F(\i;)) = (\^ + USu) T , (29) 

where u = (0,uf +1 , . . . ,uf ,0) T is a K -dimensional vector. Here the former and the latter zero 
vectors have dimensions t and K — k, respectively, and for i = t + 1, . . . , I, Uj = YlnLo Jki(^i) Uv i 
is given in Eq. 128\) . 

Proof. Under the assumption of the theorem, we have k\ = ■ ■ ■ = kt = 1 by using Lemma 16.31 
Then for any i = 1, . . . , t, 

U x S\mi) = SJ(Ux)\mi) = \S\mi), 

or in other words, S\rrii) is an eigenvector of Ux with its corresponding eigenvalue having unit 
norm. So we have P-^USlmi) = from Lemma [5lfl 2. 
On the other hand, from Eq. Q we have 

oo 

= P^ + ^PxUU'mx 

n=0 

oo 

= Px\^) + Y, P x USJ ^ J x) n S- 1 \^)x 

n=0 

oo 

= P T \4>) + P T USY,J(Ux) n v'. (30) 

n=0 

Here v' = (0, vS_j , . . . , vf) T and the zero vector has dimension t. Then the result holds by using 
Lemma 16.41 and rewriting Eq. (|30p into vector form. □ 

Although we only consider pure input states in Theorems 16.11 and 16.21 they may be used 
to calculate the outputted state F(p) of loop ([T]) for any mixed input state p by noting that 
F(p) = PiF , where p = Pil^i) {^il is the spectrum decomposition of p. 

7 Some Illustrative Examples 

To illustrate further the notions introduced and the results obtained in this paper, we consider two 
simplest classes of quantum loops. 

7.1 Single qubit loops 

Let M be an observable in the 2— dimensional Hilbert space Hi- Then we have M = mi\mi){mi\ + 
t^I^X^I) where m\,m% are the eigenvalues of M, and \rrii) is the eigenvector of M corresponding 
to mi (i = 1, 2). A single qubit loop can be written as follows: 

while (M = mi) {q := Uq}, (31) 

where U is a unitary operation on a single qubit, and i = 1, 2. Without any loss of generality we 
may assume that m\ ^ and i = 1. 
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Note that the function F defined by the loop (|3ip is from V{H2) to T>~ (span{\m2)}) ■ Since 
span{\m2)} is the one-dimensional Hilbert space, T>~ (span{\rri2)}) can be identified with the unit 
interval [0, 1]. Thus, the function F computed by the loop (|3ip is a mapping from T>(H2) into [0, 1]. 
A simple application of Theorems 14.14 15-11 and 16,21 leads to the following: 

Lemma 7.1 Let p £ D(H2) be the input state to the single qubit loop program I131\) . Then: 

1. if {m\\p\m\) = or (mi|J7|mi) = 0, then the loop \31\) terminates, and F{p) = 1; 

2. if \(mi\U\mi)\ < 1, then the loop A31\) almost terminates, and F(p) = 1; 

3. if (mi\p\mi) > and |(mi|t7|mi)| = 1, then the loop A31\) does not terminate, and F{p) = 
{m 2 \p\m 2 ). 

Now we further consider the special case that the input is a pure state. To this end, we shall 
need the following: 

Lemma 7.2 (|3J \20§ ) Each unitary operation on a single qubit can be written in the form of 
U = e ia R z ((3)R y ( , y)R z (5), where a, (3, 7 and 5 are real numbers, 

Ry(0)=( C ° S l " Sin /) 
y \ sin I cos I / 

and 

are the rotation operators about y and z axes, respectively. 

To simplify the presentation, we further suppose that the measurement is performed on the 
computational basis. Combining Lemmas 17. II and 17.21 we obtain: 

Proposition 7.1 Suppose that \ip) = ao|0) + a%\l) is the input to the single qubit loop program: 

while (q = 0) {q := Uq}, (32) 

where the loop condition (q = 0) means that the outcome of a measurement on the computational 
basis |0), |1) is 0, and the unitary operator U is given as in Lemma \7.2\ Then 

1. if ao = or 7 = (2n + l)ir for some integer n, then the loop A32\) terminates; 

2. if ao 7^ and 7 = 2mr for some integer n, then the loop \32\) does not terminate; 
3- if 7 n7T f or an V integer n, then the loop \32\) almost terminates. 

A similar conclusion holds if the guard condition (q = 0) in the loop \32\) is replaced by (q = 1), 
which means that the result 1 occurs when performing a measurement on the computational basis 
|0),|1>. 
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It is interesting to note from the above proposition that termination of the loop f|32[) depends 
only upon the parameter 7, and it is irrelevant to the other parameters a, (3 and 5. Moreover, we 
see that the loop (|32p is terminating if 7 = (2n + l)ir for some integer n, and it is (uniformly) 
almost terminating if 7 7^ nir for any integer n. 

From Lemma \7.\\ it is easy to see that in Proposition 17.11 for the cases 1 and 2, we have 
F(\ip)) = 1, and for the case 3, = |ai| 2 . 

The most frequently used single qubit gates are the four Pauli matrices: 

'=(;;)•*=(?;)• 

°0< 



Corollary 7.1 For a single qubit input state \tp) = clq\0) + ai|l), the loop ( fffffj) always terminates 
when U = X orY, it almost terminates when U = H , and it does not terminate when U = I, Z,S 
or T provided ao 7^ 0. 

7.2 Two qubit loops defined by controlled operations 

As the second example we consider a typical class of two qubit gates, namely, controlled operations. 
Suppose that U is a single qubit unitary operation. Then the controlled-?/ gate is defined by the 
following 4x4 matrix: 

where / is the 2x2 unit matrix. For a two qubit system, the measurement M on the computational 
basis 1 00), |01), 1 10) and |11) has four possible outcomes 00, 01, 10 and 11, where we use to 
indicate the measurement result q\ = i\ and 02 = ii for any i\,i2 € {0,1}. Thus, the two qubit 
quantum loop defined by controlled operation C(U) may be written as follows: 

while (M G X) { qi ,q 2 := C{U)q l ,q 2 }, (33) 

where X C {00,01,10,11}. 

The following proposition carefully examines the behavior of this loop for various choices of X 
except the trivial cases X = or X = {00, 01, 10, 11}. 

Proposition 7.2 Let pure state \tp) = aoo|00) + aoi|01) + aio|10) + an|ll) be the input of the loop 
program i33\) . Suppose that U = (Uij)jj =0 is the matrix representation of U according to the basis 
{|0), |1)}, that is, Uij = (i\U\j) for any i,j E {0, 1}. 



the Hadamard gate: 

H 



111 



V2\l -1 
the phase gate: 

and the f gate: 

' 1 

expCf) 

Applying Proposition 17. II to these gates, we obtain: 
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1. IfX = {00}, thenpNT = \aoo\ 2 , F(\ip)) = a i|01) + ai |10) + an|ll), the loop fiJ3\) terminates 
if aoo = 0, and it does not terminates if aoo ^ 0. 

2. IfX = {01}, thenp N T = \a i\ 2 , F(\ip)) = o o|00)+aio|10) + an|ll), the loop {HP terminates 
if a oi = 0, and it does not terminates if aoi ^ 0. 

3. Let X = {10}. If a,io = or Uqq = 0, then the loop fiJty) terminates. If aio = or \Uoo\ < !• 
then it almost terminates, and 



If a io 7^ and Uqo = 1, then it does not terminate, and F(\tp}) = aoo|00) +aoi|01) + an|ll). 

4- Let X = {11}. If an = or Un = 0, then the loop fifty) terminates. If an = or \Un\ < 1, 
then it almost terminates, and 



If a n 7^ and Un = 1, then it does not terminate, and F(\ip)) = aoo|00) +aoi|01) + aio|10). 

5. If X = {00,01}, then p NT = |a 00 | 2 + |a i| 2 , = ai |10) + a n \ll), the loop [Sty) 
terminates if aoo = aoi = 0; an d it does not terminate if aoo / or aoi ^ 0. 

6. If X = {10,11}, then p NT = \a 10 \ 2 + |an| 2 , F(\ip)) = a 00 |00) + a O i|01), the loop fiJty) 
terminates if aio = o-n = 0, and it does not terminate if aio ^ or an 7^ 0. 

7. Let X = {00, 10}. Then we have: 



F(\ip)) £ V (span{\01), |11}}) is given as follows: for the case of \Uoo\ = 1, F(\ip}) = 
aoi|01) +an|ll), and for the case of \Uqq\ < 1, 



If a oo = 0, and aio = or Uoo = 0, then the loop fifty) terminates, if ooo = 0, and aio = or 
| fool < 1; then it almost terminates, and if aoo 0; or °io 7^ and \Uqq\ = 1, then it does 
not terminate. 

8. Let X = {00, 11}. Then we have: 




£D(span{|00),|01),|ll)}). 




£D(span{|00),|01),|10)}). 






23 



F(\ip)) G V (span{\01), |10)}) is given as follows: for the case of \Un\ = 1, F(\ip}) 
aoi|01) +aio|10), and for the case of\Uu\ < 1, 

w» = ( l aol S , a TfZ p 

V °io a oi l a iol + l a n| 

If aoo = 0, and an = or U\\ = 0, the the loop t33\) terminates, if aoo = and \Uu\ < 1, 
or aoo = and an = 0, then it almost terminates, and if aoo 7^ 0> or a n 7^ anc ^ l^nl = 1; 
then it does not terminate. 



9. Let X = {01, 10}. Then we have: 



Pnt 



\aoi\ 2 , if\U 00 \ < 1, 

Ki| 2 + |aio| 2 , if \U o\ = 1. 



€ 2? {span{ | 00), 1 11)}) is given as follows: for the case of \Uoo\ = 1, -F(IV')) 
aoo|00) +an|ll), and for the case of\Uoo\ < 1? 

V a n a oo l a iol + l a n| 

If a oi = 0; an d a w = or Uoo = 0, the the loop [&3]) terminates, if aoi = 0, and |f/oo| < 1 or 
a io = 0, then it almost terminates, and if aoi 7^ 0, or aio 7^ and |E/oo| = 1; i/ien i£ does not 
terminate. 



10. Let X = {01, 11}. Then we have: 



PNT 



Ki| 2 , if\Uu\ < 1, 

Ki| 2 + |an| 2 , if\Uu\ = 1. 



F(|^)) € £> (span{ 1 00), 1 10)}) is given as follows: for the case of \Un\ = 1, F(\tJ))) 
ooo|00) + aio|10), and for the case of\U\\\ < 1, 

FW) = ( , a Tfla P 

V a 100oo l a iol + l a ll| 



If a oi = 0, and an = or Uu = 0, the the loop \33\) terminates, i/aoi = 0, and \U\\\ < 1 or 
an = 0, then it almost terminates, and if aoi 7^ 0, or an 7^ and \U\\\ = 1, then it does not 
terminate. 



11. Let X = {00,01, 10}. Then we have: 

PNT = 

and F(\ijj)) S T>~ (span{\ll)}) = [0, 1] is given by 

F(H)) 



aool 2 + Ki| 2 , if \Uoo\ < 1, 

oool 2 + Ki| 2 + Ko| 2 , if\Uoo\ = l, 



aiol 2 + |an| 2 , if \U o\ < 1, 
on| 2 , if IE/00 1 = 1- 



If a oo = 001 = 0, and aio = or [Too = 0, then the loop [33]) terminates, if aoo = &oi = 0, 
and \Uoo\ < 1 or aio = 0; then it almost terminates, and if aoo 7^ 0, or aoi 7^ 0, or aio 7^ 
and \ Uoo\ = 1, i/ien i£ does not terminate. 
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12. Let X = {00,01, 11}. Then we have: 

(\a 00 \ 2 + |a i| 2 , if\Un\ < 1, 

^aoor + laoil+lanl, *fpii\ = 1, 

and F(|V>)) G £>~(spare{|10)}) = [0,1] is given by 



|aio| 2 + |an| 2 , if\U n \ < 1, 
|aio| 2 , if \Uu\ = 1. 



If a oo = ooi = 0, and an = or Un = 0, £/ten the loop {33\) terminates, if aoo = «oi = 0, 
and \U\i\ < 1 or an = 0, then it almost terminates, and if aoo 7^ 0, or oq% ^ 0, or an / 
and \ Un\ = 1; then it does not terminate. 

13. LetX = {00,10,11}. Thenp NT = |a o| 2 +|ai | 2 +|an| 2 andF(\i/))) = \a i\ 2 G V~ (span{|01)}) 
[0, 1]. If aoo = a io = a n = 0, then the loop [33\) terminates, otherwise it does not terminate. 

14. LetX = {01,10,11}. Thenp NT = |a i| 2 +|ai | 2 +|an| 2 andF(\i/))) = |a 00 | 2 G V~ (span{\00)}) 
[0, 1]. If aoi = a io = a n = 0, i/ien the loop [33\) terminates, otherwise it does not terminate. 

Note that termination of the loop (|33p is irrelevant to the unitary operator U, and it only 
depends on the input state |V>> when X = {00}, {01}, {00, 01}, {10, 11}, {00, 10, 11} or {01, 10, 11}. 
For the other cases, termination of the loop defined by the CNOT gate is summarized in the 
following: 

Corollary 7.2 Suppose that C(U) is the CNOT gate C(X), where X = NOT is the second Pauli 
gate. 

1. Let X = {10} or {11}. Then the loop \33\) always terminates. 

2. Let X = {00, 10} or {00, 11}. Then the loop 133\) terminates if aoo = 0, otherwise it does not 
terminate. 

3. Let X = {01, 10} or {01, 11}. Then the loop 133\) terminates if aoi = 0, otherwise it does not 
terminate. 

4- LetX = {00,01, 10} or {00, 01, 11}. Then the loop i33\) terminates if aoo = a oi = 0, otherwise 
it does not terminate. 



8 Conclusion 



Exploitation of the full power of loop construct in quantum computation requires a deep under- 
standing of the computational mechanism of quantum loop programs. In this paper, we introduced 
a general scheme of quantum loop programs, the behaviors of quantum loops are carefully analyzed, 
including termination, almost termination, and sensitivity, and a matrix-summation representation 
of the function computed by a quantum loop is found. 

This paper is merely an initial step toward a thorough understanding of quantum loop programs, 
and many important problems remain open. First, the bodies of quantum loops that we considered 
are unitary transformations. If a quantum loop is allowed to be embedded into another quantum 
loop, then as was observed in Section 3, the body of the latter loop is not a unitary operator 
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but a super-operator in general. Therefore, it is an interesting topic for further studies to find 
conditions for termination and almost termination of quantum loops in which the loop bodies are 
super-operators. Second, we demonstrated the expressive power of quantum loops by presenting 
a loop description of quantum walks. It would be very interesting to find more computational 
problems that cannot be expressed or solved without quantum loops. In general, the study of loop 
programs is a very important area of computer programming methodology. Reconsideration of 
some fundamental problems from this area, say loop invariants and proof rules, in the quantum 
setting is a great challenge. 
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